Butterfly Diagram and Activity Cycles in HR 1099 
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ABSTRACT 

We analyze photometric data of the active RS CVn-type star HR 1099 for the years 1975- 
2006 with an inversion technique and reveal the nature of two activity cycles of 15-16 yr and 
5.3±0.1 yr duration. The 16 yr cycle is related to variations of the total spot area and is coupled 
with the differential rotation, while the 5.3 yr cycle is caused by the symmetric redistribution of 
the spotted area between the opposite stellar hemispheres (flip- flop cycle). We recover long-lived 
active regions comprising two active longitudes that migrate in the orbital reference frame with 
a variable rate because of the differential rotation along with changes in the mean spot latitudes. 
The migration pattern is periodic with the 16 yr cycle. Combining the longitudinal migration 
of the active regions with a previously measured differential rotation law, we recover the first 
stellar butterfly diagram without an assumption about spot shapes. We find that mean latitudes 
of active regions at opposite longitudes change antisymmetrically in the course of the 16 yr cycle: 
while one active region migrates to the pole, the other approaches the equator. This suggests a 
precession of the global magnetic field with respect to the stellar rotational axis. 

Subject headings: stars: activity, imaging, rotation, spots, individual (HR 1099) — binaries: 
spectroscopic — techniques : photometric 



1. Introduction 

HR 1099 (= V711 Tau) is the most exten- 
sively studied RS CVn system, consisting of a 
Kl IV primary and a G5 IV secondary in a 2.8 
day orbit. The primary, being the faster rota- 
tor, is the first cool active star whose surface was 
mapped with the Dopp ler imaging (DI) technique 
([Vogt fe Penrodl Il983l ) . It has a record 46 sur- 
face maps obtained between 1981 and the present 



area of the polar spot and its extensions may vary 
cyclically and correlate with th e mean brightness 
of the system (jVogt et al.lll999h . 



( Donati et al.lll992l.[2003l:lDonat 



20031: iPetit et al 
star feature 



2004) 



1999HStrassmeier fc BartusfcOOC ; 



1999; Vogt et al 



Garcia- Alvarez et al 



Almost all images of the 
a prominent, cool, polar cap along 
with a variety of high- and low-latitude spots. The 



The photometric variability of HR 109 9 was dis- 
covered in 1975 by lLandis fc Hall ( 1976?) and has 
been followed regularly for more than 30 years. 
The dat a before 199 3 were collected and ana- 
lyzed by iHenrv et al.1 (|l995l ). They found peri- 
odicities of 5.5±0.3 and 16±1 yr, respectively, in 
the star's mean and maximum brightness, suggest- 
ing the presence of activity cycles. Using a model 
with two circular spots, they identified 15 differ- 
ent dark spots within the 18 yr span of the data 
with lifetimes averaging between 1.5 and 2.0 yr. 



1 



The spot rotation periods ranged from 2.8245 to 
2.8412 days and bracketed the orbital period by 
—0.47% to +0.12 % days, implying significant dif- 
ferential rotation. Vogt et all (|l999l ) measured the 
differential rotation from DI and found that it is 
of opposite sign and about a factor of 56 less than 
for the Sun . With newer stellar images, however, 
Petit et al. I (|2004h argued that the differential ro- 
tation law for HR 1099 is of the same sign as the 
Sun and a factor of 2 larger than determined by 
Vogt et al. ( 1999h . The ran ge of polar and equ ato- 
rial periods determined bv lPetit et al. (200J) en- 



comp asses all periods de rived from the photome- 



try of iHenrv et all (|l995l ) 



Mohin fc Raveendranl ( 1993f ) interpreted the 



photometric variability of HR 1099 as the result 
of two active longitude regions migrating in the or- 
bital reference frame with constant but very differ- 
ent rates. Two active longitudes (ALs) separated 
by about 180° and locked in synchron ism with the 
orbita l period have been surmised bv lHenrv et al. 
(|l995l ). In comparison with other RS CVn-type 



stars, however, the ALs on HR 1099 appeared very 
broad, 110°-140°, implying smearing by differen- 
tial rotation. 

In the present Letter we analyze all photo- 
metric data collected up to date for HR 1099, 
previously published and newly obtained. Em- 
ploying a technique based on inversions of stel- 
lar l ight curves into images of the spot filling fac- 
tor (|Berdvugina et al. I l2002h . we recover locations 
of active regions and search for long-lived ALs. 
Taking into account the effect of differential rota- 
tion, we are able to recover a stellar butterfly dia- 
gram, which, in con trast to previous studies (e.g., 
Katsova et alT2 0031 , does not involve assumptions 
about spot distributions, shapes, and numbers. 

2. Observations 

We have analyzed the Johnson V photometric 
observations collected earlier from the literature 
and from the Tl 0.25 m and T3 0.40 m automatic 
photo metric telescopes (APTs) by IHenrv et~al 



(|1995I ). that span the years 1975-1993. We have 
also included new observations for the years 1994- 
2006, acquired with the T3 APT. The Tl and T3 
APTs, along with several others, are housed at 
Fairborn Observatory in southern Arizona. Fur- 
ther details on the operation of the APTs and 



the red uction of the data can be found in IHenrv 
(|l995al fbl) and lEaton et all (|2003l ). The complete 
set of Johnson V differential magnitudes is listed 
in Table [1] and plotted in Fig. [TJ The differential 
magnitudes have been converted to apparent mag- 
nitudes in Fig. [TJj with a V magnitude of 4.28 for 
the comparison star 10 Tau (= HR 1101). The 
faint, close, visual companion of HR 1099 was al- 
ways included in the photometer diaphragm and 
so is included in the light of HR 1099. 

The V observations were divided into 109 sub- 
sets within which the light curve was observed to 
be stable over each time interval. The average 
length of each subinterval was <40 days. Julian 
Dates of the observations were converted to orbital 
phases with the ephemeris 



conj 



= 2, 442, 766.08 + 2.83774.B 



(1) 



adopted from iFekell ()1983f ). where zero phase cor- 
responds to the conjunction with the Kl star clos- 
est to the Earth. The error in the period of 
10~ 5 days results in a phase uncertainty of only 
0.04 af ter 30 years. Accord ing to the Fourier anal- 
ysis of I Henry et al.l ( 19951 ). the star exhibits the 
presence of a small ellipticity effect with a full am- 
plitude of 0.016 mag in V . This effect was removed 
from each data set before the light-curve inversions 
were performed. 

3. Activity cycles and active longitudes 

For each of the 109 light curves, we measure 
the maximum, mean, and minimum V magnitudes 
as well as the peak-to-peak amplitude in V. Pe- 
riod anal ysis of these values confirms the cycles 
found bv iHenrv et all (Il995f ). The mean bright- 
ness variations show a cycle of ~15 years (Fig.UJi), 
and similar cycle lengths of 15-16 years are seen 
in the maximum and minimum magnitudes. The 
peak-to-peak amplitude in V varies with a cycle 
of 5.3±0.1 yr (Fig. Q})). This cycle also has an in- 
fluence on variations in the maximum, mean, and 
minimum V magnitudes. Similarly, the 16 yr cy- 
cle also influences the amplitude variations, as the 
largest amplitude of 0.2 mag was observed at the 
maxima of the spottedness in 1979 and 1994. 

To determine the large-scale spot distribution 
on HR 1099, we apply an inversion technique to 
the light curves using a two-tem perature approx- 
imation (jBerdvugina et al.ll2002i ). The inversion 
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of the light curve results in the distribution of 
the spot filling factor over the stellar surface, i.e., 
a stellar image (Fig. [2|). We adopted values of 
4700 and 3500 K for the photosph ere and spot 



tempe ratures, respectively, following lHatzes et al 
(Il996h . Also, we assumed that the star is spot- 
less at V = 5.6 mag, which is about 0.05 mag 
brighter than the brightest magnitude observed for 
HR 1099. The inclination of the rotational axis to 
the line of sight is chosen to be 40°; inclinations 
in the range of 33°-45° have been used for the 
previous DI mentioned above. 

From the stellar images, we recover spot lon- 
gitudes with the maximum sp ot filling factor 



see details in iBerdvugina et al.l 120021 ) . The ac 



curacy of the spot longitudes depends on the 
broadness of the light curve minimum and is 
on average 0.05 in phase. When two minima 
are seen in the light curves, two spot concen- 
trations are recovered in the images. The time 
dependence of the spot phases is shown in Fig. [T^ 
(the data are available on request). Our derived 
spot phases constitute two sequences that repre- 
sent two migrating active longitude regions sepa- 
rated on average by 180°, similar to the previous 
studies of HR 1099 iMohin fc Raveendranlll993t 
Henry et al. 19951 : lLanza et al.ll2006f) and to other 



RS CVn-type stars (IBerdvugina fc T uomincn 



1991 iFekel fc Hemvll2005h . lLanza et all (|2006h 
related the long-term trend in the migration of 
the ALs with binary orbital period variations. 

We notice that the two active regions regularly 
approach each other and merge into a large sin- 
gle region. Subsequently, the single region divides 
back into two spot regions. Thus, we are seeing the 
crossing of the migration paths of the two ALs. 
In the presence of differential rotation, such be- 
havior suggests variations in mean spot latitudes. 
Three intersections were clearly observed in 1979, 
1990, and 1995. At the same times, the amplitude 
of the brightness variations reached a maximum 
(Fig. [TJ>) , which implies a rearrangement of spots 
onto one stellar hemisphere and supports the idea 
of merging active regions and crossing ALs. The 
intersections in 19 79 and 1990 a re also seen in the 
reconstruction bv lLanza et al. ( 2006f ). However, 
they missed the one in 1995 due to lack of data. 
The shorter (~5 yr) and longer (~11 yr) time in- 
tervals between the intersections comprise the 16 
yr cycle. 



In Fig. [Tb, the two migration paths are plotted 
with two fits to the observed phases. The fits are 
obtained with a two-harmonic Fourier series and 
a long-term trend, which together are able to re- 
produce the nonsymmetric shapes of the curves. 
The assignment of the active regions to one or the 
other AL has been done under the assumption that 
the active region continues its migration with the 
same rate after the intersection of the ALs, which 
is natural in the presence of differential rotation. 
Although the fits have been calculated separately 
for the two ALs, both reveal similar shapes that 
are mirrored with respect to each other in phase. 
The migration patterns appear to be periodic with 
a 16±3 yr cycle. The standard deviation of the fits 
is as small as 0.06 phase units. An extrapolation 
of the fits suggests that the next crossing of the 
active longitudes will occur in 2007. Observation 
of the ALs crossing in 2007 with an increase in 
photometric amplitude will be a good test for the 
present model. 

In addition to the migration patterns and 
their intersections, abrupt switches of the dom- 
inant activity between the ALs occurred in 1988, 
1992-1993, 1998-1999, and 2005, all at the times 
when the amplitude of the light curve was near 
the minimum. Such events periodically occur in 
RS CVn- and FK Corn-type stars and are known 
as flip-flop cycles (jBerdvugina &: Tuominenl 1998t 
iKorhonen etaD |2002j). In HR 1099, the flip-flops 
occurred when the active regions at opposite lon- 
gitudes became about the same size; that reduced 
the amplitude of the light curve, and the dominate 
activity switched to the opposite longitude. Thus, 
flip-flops in HR 1099 are apparently governed by 
the 5.3 yr cycle but appear to be disturbed by the 
differential rotation and intersections of the active 
longitudes. 

4. The butterfly diagram 

The shape of the spot migration curves for 



HR 1 099 is reminiscent of the solar ALs (jBerdvugina fc Usoskin 
2003). On the Sun, such a migration pattern is 
determined by changes of the mean spot latitudes 
along with the surface differential rotation. As- 
suming that a similar relation holds for HR 1099, 
we can recover mean spot latitudes on its sur- 
face, i.e., create butterfly diagrams for the active 
longitudes in the stellar hemisphere visible above 
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the stellar equator. For this purpose we use the 
differential rotation law: 



fi = O eq - Afi sin 2 i(>, 



(2) 



with the parameters estimated by IPetit et"aTI 
(|2004l ) from DI of HR 1099 for the year 2002, hav- 
ing the smallest variances: fi cq =2.2241±0.0004radday _1 
and Afi=15.2±0.8mradday _1 . By observing a 
phase shift of the ALs during a given time inter- 
val, we can determine an instant angular velocity 
fi, and, using Eq. ([2]), we can find the correspond- 
ing latitude ip. To increase accuracy, we determine 
phase shifts from the harmonic fits to the migra- 
tion paths. The results are shown as solid and 
dashed lines in Fig. [3jz & b. As expected, mean 
spot latitudes vary with the 16 yr cycle for both 
active longitudes. 

The extent and shape of the butterfly dia- 
gram recovered from the AL migration can be 
compared with the measurements of spot lati- 
tudes from DI. From the 46 published images, we 
measured latitude extents of the recovered spots 
and calculated weighted mean latitudes for each 
of the two opposite ALs. These are shown in 
Fig. [3Jz & b. The mean DI latitudes and the 
curves obtained from the AL migration are in 
good agreement, thus supporting our interpreta- 
tion of the photometry. To estimate a possible 
range of the parameters in Eq. ([2]), we searched 
for the best fit (via x 2 minimization) to the mean 
DI latitudes combined with the angular veloc- 
ity fi measured from the AL migration. The re- 
sult is shown in Fig. 2J It reveals a very elon- 
gat ed error elli p soid e ncompassing measurements 
by IPetit et all ((200l for different years. Our 
best-fit parameters, fi cq =2.230±0.003radday _1 
and Afi=22±4mradday -1 , indicate a somewhat 
stronger differential rotation on HR 1099 than the 
best parameters bv IPetit et al. ( 2004 ). although 
the two measurements appear to be statistically 
equivalent. 

The shape of the butterfly diagram for HR 1099 
is reminiscent of the Sun's. The direction of the 
latitude migration, however, is different for the 
two ALs. For the active region shown in Fig. QJ 
with a solid line, preferred spot appearances first 
move relatively rapidly (for 5 years) to higher lat- 
itudes, up to 70°, and then slowly (for 11 years) 
move back to mid latitudes, down to 40° (Fig.[3}i). 
For the other active region (shown in Fig. with a 



dashed line), the behavior of the spot latitudes is 
just the opposite. The latitude shifts are there- 
fore always opposite for the two ALs. During 
1982-1987, both ALs were at about the same lat- 
itude, and the migration rate remained constant. 
In 1988-1993, one active longitude was occupied 
by large high-latitude (polar) spots, while low- 
latitude spots appeared around the opposite active 
longitude. After 1997, the dominant active regions 
again moved to about the same latitude, thus com- 
pleting the 16 yr cycle. The opposite latitudinal 
migration of the two long-lived regions could be 
responsible for the contradi ctory measuremen ts of 
the differential ro tation bv lVogt et al. ( 19991 ) and 
Petit et all (|2004 . 



Discussion 



As noted bv iHenrv et al. ( 1995 ). if the 16 yr 
cycle is akin to the 11 yr sunspot cycle, one 
might anticipate other signatures of such a cy- 
cle. Indeed, the 16 and 5.3 yr activity cycles 
in HR 1099 are s imilar to the sunspot cycles 
of 11 and 3.7 yr (|Berdvugina fc Usoskinl [2003? ) 
and the cycles observed on y oung solar analogs 
([Berdvugina fc Jarvinenll2005r ). The longer cycle 
in HR 1099 reflects the variation of the total spot- 
tedness and is coupled with the differential rota- 
tion. This cycle may be associated with an oscil- 
lating axisymmetric dynamo mode in which global 
polarity reversals occur in cycles of 22 years on the 
Sun or 32 years on HR 1099. No reversals have yet 
been observed on HR 1099 des pite 12 years of spec - 
tropolarimetric observations (jDonati et al. 2003). 
Comparing the butterfly diagrams of the Sun and 
HR 1099, we suppose that a global polarity rever- 
sal on HR 1099 could have happened around 1990, 
just before the regular spectropolarimetric obser- 
vations by Donati et al. began. The next one can 
be expected in 2007. 

The shorter activity cycle on HR 1099 is related 
to the flip-flop phenomenon; i.e., it reflects sym- 
metrical redistribution of the spotted area between 
opposite stellar hemispheres. In the light curves, 
the flip-flops are seen as a strengthening of the sec- 
ondary minima and a diminishing of the primary 
on e. In the Dopp ler im ages of HR 1099 obt ained 
by iDonatil (|l999h and iDonati et all (|2003l ). the 
flip-flop in 1992-1993 is seen as an abrupt change 
of the phase towards which the polar spot is tilted 
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by about half a rotational period, while the flip- 
flop in 1998-1999 appears as a progressive shift 
of the polar spot extension to the opposite longi- 
tude. The flip-flop cycles can be associated with a 
non-axi symmetric dynamo mode, which gives rise 
to ALs (jFluri fe Berdvuginall2004l) . Similar to the 
Sun and young solar analogs, the flip-flop cycle on 
HR 1099 is about one-third of the length of the cy- 
cle regulating the total spottedness, which implies 
a certain relation between the non- and axisym- 
metric dynamo modes. 

The fact that the butterfly diagrams for the two 
ALs are antisymmetric, i.e., that latitude shifts 
are always opposite in the opposite stellar hemi- 
spheres, implies that the axis perpendicular to 
the line connecting the ALs changes its direction 
with respect to the rotational axis by about ±15° 
during the 16 yr cycle. Thus, the two axes are 
separated the most in 1990-1991 and 1995-1996 
and aligned in 1982-1987 and 1998-2003, sug- 
gesting the precession of the nonaxisymmetric dy- 
namo mode during the 16 yr cycle. Such a pre- 
cession of the magnetic field was also suspected 
in the Ze eman-Doppler imag es of HR 1099 ob- 



tained by iDonati et alj (|2003f ) , who noticed that 
the rings of the azimuthal field got progressively 
better centered in 2000-2002, while the fraction 
of the magnetic energy attributed to the merid- 
ional field component was regularly increasing in 
1992-1997 and decreasing starting from 1998 up to 
now. They assumed that such behavior suggests 
that the main axis of symmetry of the azimuthal 
field component is varying in orientation with time 
on the timescale of at least one decade, which is 
in agreement with our findings. 

Our results suggest that a few interesting events 
in the magnetic activity of HR 1099 may happen 
in 2007. For instance, an extrapolation of the fits 
to the ALs predicts their next crossing to happen 
during this time. Also, we can expect the mean 
spot latitude of the dominant active region to be 
different again and the rings of the azimuthal field 
to be decentered. In addition, a global polarity 
reversal may finally be observed. The detection of 
such events would be of great importance for un- 
derstanding stellar magnetic cycles and their un- 
derlying dynamo processes. 

We thank Lou Boyd at Fairborn for making 
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search at Tennessee State University was sup- 
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1975. 19B0. 1985. 1990. 1995. 2000. 2005. 

Year 

Fig. 1— Activity patterns on HR 1099. (a) All 
V-magnitude observations an alyzed in this Letter . 
The data before 1993 are from Henry et al. ( 1995h . 
Cyclic variations of the mean brightness are re- 
vealed with a harmonic fit (solid curve) with a 
period of ~15 years, (b) Peak-to-peak amplitude 
of the light curves (/illed circles) . A harmonic fit 
(solid curve) reveals the period of 5.3±0.1 yr. (c) 
Phases of the spots recovered with the light-curve 
inversions. Filled circles denote the larger or sin- 
gle active region observed at a given time, while 
open circles indicate a second, smaller region if ob- 
served. Migration paths of the ALs are indicated 
with solid and dashed curves, which are harmonic 
fits to the measurements with a long-term trend. 
The migration pattern is apparently repeated in 
about 16 years. 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Table 1 
HR 1099 Photometric Data. 



HJD 


AV 


2,442,719.8650 


1.600 


2,442,720.9300 


1.510 


2,442,721.9050 


1.570 



* Table 1 is published 
in its entirety in the elec- 
tronic edition of the As- 
trophysical Journal. A 
portion is shown here 
for guidance regarding its 
form and content. 




Phase Phase 

Fig. 2. — Examples of the spot filling factor distri- 
bution on HR 1099 (Zeft) and corresponding light 
curves (right). The stellar surface was modeled 
with a grid of 10°xl0°. The spot filling factor 
at each grid point indicates relative contributions 
from two temperature components, the hot pho- 
tosphere and cool spots. In the images, the darker 
area implies a higher spot filling factor. Observed 
light curves are shown as crosses and the fit as the 
solid line. In contrast to direct modeling of light 
curves, the inversions do not involve any assump- 
tions about spot shapes, numbers, or latitudes. 
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Fig. 3.— Butterfly diagram of HR 1099. (a) 
Mean latitudes recovered from the fit to the ac- 
tive longitude shown in Fig. [TJ; as the solid curve. 
(6) Mean latitudes recovered from the fit to the 
active longitude shown in Fig. QJ: as the dashed 
curve. In both plots, the mean latitudes of spots 
in the opposite hemispheres measured from DI are 
shown with different symb ols corresponding to th e 



follow ing sources: VP 83 dVogt fc Penrod 
V99 (IVoet et alJll999h . D9 2 dDonati et al ' 



1983), 



1992) 



D99 (lDonatil ll999l). SB00 dStrassmeier fe Bartusl 



2000). G03 (IGarria- Alvarez et alj 120031) . D03 
(jDonati et aT] l2003h . P04 (|Petit et al.ll2004h . and 
P06 (P. Petit 2006, private communication). Stan- 
dard deviations of the mean DI latitudes (if sev- 
eral spots were observed) are shown as vertical 
dashes. The thick curves are the solutions using 
th e best parameters of the differential rotation law 
bv lPetit et ail (2004). The thin curves indicate the 
±1(7 interval for the best parameters obtained in 
this Letter (filled circle in Fig. [4]) . The two solu- 
tions are statistically equivalent. 




eq , rad/day 

Fig. 4. — The \ 2 surface contours for the two pa- 
rameters in Eq. (|2|) calculated from the AL migra- 
tion and the DI mean spot latitudes. The x 2= l 
level is emphasized by a thick contour. The best- 
fit solution is ind icated by the fille d circle. The 
measurements by Petit et al. ( 20041 ) are shown as 
open circles. The strongly elongated shape of the 
minimum implies that the two parameters are al- 
most linearly dependent. The distribution of the 
measurements along the major axis of the error el- 
lipsoid suggests that they reflect errors in DI spot 
latitudes rather than true variations of the differ- 
ential rotation. 
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